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Fermi-surface deformations in strongly correlated metals, in comparison to results from 
band-structure calculations, are investigated. We show that correlation-induced interband 
charge transfers in multi-orbital systems may give rise to substantial modifications of the 
actual Fermi surface. Depending in particular on the relative strength of the crystal-field 
splitting and of the Hund's exchange coupling, correlations may either reinforce orbital po- 
larization or tend to compensate differences in orbital occupancies, as demonstrated by 
investigating a 2-band Hubbard model in the framework of dynamical mean field theory 
(DMFT). The physical implications of such interorbital charge transfers are then explored 
in two case studies: BaVS3 and Na^CoC^. By means of the DMFT in combination with the 
local density approximation (LDA) to density functional theory (DFT), new insights in the 
underlying mechanism of the metal-to-insulator transition (MIT) of BaVS3 are obtained. 
A strong charge redistribution in comparison to LDA calculations, i.e., a depletion of the 
broader A\ g band in favor of the narrower E g bands just above the MIT is found. In ad- 
dition, the intriguing problem of determining the Fermi surface in the strongly correlated 
cobaltate system Na^OoC^ is discussed. 

§1. Introduction 

In recent times, theoretical investigations of systems with strong electronic corre- 
lations have become feasible in realistic settings by using the material-specific band 
structure input as a single-particle reference systemP This also raises new ques- 
tions, in particular concerning effects stemming from the multi-orbital character of 
realistic materials. As most physical properties of a specific material in the metallic 
state are determined by the electronic states close to the Fermi level, one of the 
key questions is how these states are affected by strong Coulomb correlations. It is 
a well-documented fact, e.g., in cuprates, that the competition of localization and 
itinerancy can result in the appearance of distinct areas on the Fermi surface with 
strikingly different physical properties of the low-energy excitations. In this work, 
we concentrate on a more basic aspect, namely how correlations may lead to a mod- 
ification of the shape of the Fermi surface itself, due to interorbital charge transfers. 
These modifications may result in a major change of the low-energy physics, e.g., by 
inducing new instabilities. 

The track record of density functional theory (DFT) in local (spin) density 
or generalized gradient approximation, L(S)DA or GGA*\ promoted the effective 
single-particle approach to a standard tool for electronic structure investigations. 
Despite the many successes, this approach is generally inadequate for strongly cor- 
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related systems, and qualitative as well as quantitative discrepancies are well doc- 
umented. However, surprisingly, for many correlated metals, such as S^RuC^P 
optimally doped cuprates, and several heavy- fermion systems,'^ the Fermi surface 
predicted by LDA is apparently in rather good agreement with experimental data. 
Of course, there are also famous counter examples, such as for instance the LDA- 
predicted hole pockets in fcc-Ni which are absent in experiments^ Nevertheless, 
the overall encouraging results for the Fermi surface obtained by LDA calculations 
provided 'empirical' confidence that the Fermi surface is little affected by strong 
correlations in many materials. In this paper, we want to draw attention to the 
fact that there are indeed cases where electronic correlations do have a significant 
influence on the relative orbital occupancies, and on the Fermi-surface shape. This 
may result in a dramatic change of the low-energy physics, e.g., by allowing for new 
nesting possibilities. 

Consider first the simple case of a perfectly cubic material involving a degen- 
erate t2 g multiplet. The effective single-particle (Kohn-Sham) hamiltonian can be 
diagonalized to yield the bands e^ - with orbital index m. Assuming a local (i.e. 
k- independent) self-energy matrix, cubic symmetry implies: E mm i{oS)=8 mm iE{oS). 
Hence, the Fermi surface of the interacting system corresponds to the fc-points such 
that: s^=^ — ^(0) (with \i the chemical potential), implying that all Fermi-surface 
sheets are shifted in the same manner by correlations. However, Luttinger's theo- 
rem® implies that the total volume of the Fermi surface is unchanged by interac- 
tions, and coincides with the total number of electrons. As a result, in this simplest 
case, each Fermi surface sheet cannot be changed by correlations and the relation 
[i — U(0)=[io(n) must be satisfied (with fiQ the chemical potential of the reference 
system corresponding to the given total electron number). This teaches us that, 
for perfect cubic symmetry and a degenerate multiplet, the Fermi surface can only 
be changed by correlations if the self-energy has strong k-dependence. In contrast, 
when orbital degeneracies are lifted due to strong anisotropies or crystal-field effects, 
even a local self-energy matrix can in principle induce Fermi-surface deformations 
due to its multi-orbital (non-scalar) structure. The interacting Fermi surface is now 
determined by: det[//<5 mm / — H^ m — I7 mm /(0)]=0. Even if the off-diagonal compo- 
nents of the self-energy are negligible, each Fermi-surface sheet can be affected in 
a different manner according to: e^=/U — S mm (0). Note that Luttinger's theorem 
does not apply to each sheet separately, i.e., the partial orbital occupancies n m do 
not necessarily correspond to the respective volume of the different sheets (even for a 
diagonal £). However, physical intuition suggests that large correlation- induced in- 
terorbital charge transfers (i.e., changes of the partial n m for a fixed n=J2 m n m), will 
be associated with strong Fermi-surface deformations, as confirmed by the studies 
below. 

In this paper, we first explore the relevant processes and mechanisms that de- 
termine interband charge transfers in multi-orbital systems. Such charge transfers 
have been recently discussed by several authors 0J3J3J3i3Z3Jn} With the vanadium 
sulfide compound BaVS3 and the puzzling sodium cobaltates, we then discuss two 
examples of systems where such repopulations may drive Fermi-surface deformations 
that play an essential role for the relevant physical properties. 
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§2. Multi-orbital effects in correlated systems 

In real materials, the particularities of the band structure can have substantial 
influence on the effect of correlations. Crystal- field splittings and anisotropies giving 
rise to partial densities of states with different bandwidths and shapes lead in general 
to different renormalizations of the different bands by the Coulomb interactions. The 
net results are charge transfers between different subbands induced by the Coulomb 
interactions and a modified Fermi surface as compared to an effective single-particle 
reference system, e.g., as given by the LDA. The shape of the Fermi surface and the 
different subband fillings are in general determined by the interplay of crystal-field 
splitting, bandwidth ratios (or more generally the shapes of the partial densities of 
states), hybridizations, Coulomb interactions - in particular Hund's rule coupling - 
as well as the global filling of the system. On general grounds, the electrons tend 
to either avoid the high cost of Coulomb interactions or to compensate it by enlarg- 
ing the kinetic energy. Since the Coulomb interactions are smaller for electrons in 
different orbitals, uniform occupations of the subbands are naturally favored by the 
multi-orbital interaction vertex. The difference between the intraorbital Coulomb 
interaction U and the interorbital one U' is mathematically related to the Hund's 
rule integral J that describes the lowering of the Coulomb interaction by exchange 
between electrons with equal spins. *' Thus, a strong Hund's rule coupling will 
generally lead to rather uniform partial occupations of the different subbands. The 
crystal-field splitting A, on the other hand, tends to orbitally polarize the system. 
As shown by Manini et a/. 131 in the framework of a two-band Hubbard model, the 
relevant energy scale that has to be compared to the crystal-field splitting is the 
renormalized bandwidth, not the bare one. In a system close to the Mott transition, 
small crystal-field splittings can therefore be sufficient to induce a strong orbital 
polarization. Recently, this effect was showrP to be crucial to understand the insu- 
lating oxides LaTiOs and YTiOs. 

The overall orbital charge distribution is thus the result of a competition be- 
tween Hund's coupling which tends to level out differences in the orbital occupa- 
tions and crystal-field splittings which favor orbital polarization. This interplay 
is influenced by the band structure of the material: in narrow band systems or 
for large Coulomb interactions U the crystal-field splittings are likely to dominate 
and the system may eventually be driven into an orbitally polarized state. More- 
over, in strongly anisotropic systems involving bands with very different bandwidths, 
correlation- induced band narrowing can be different in the various bands, also af- 
fecting the magnitude of the charge transfer. The final state is thus a subtle balance 
between the Hund's rule coupling, the effective crystal-field splitting and the ef- 
fective bandwidths of the subbands, the latter two renormalized by the Coulomb 
interactions. 

To substantiate these qualitative expectations, we have studied the interplay of 
the crystal-field splitting and Hund's rule coupling as a function of U in an anisotropic 
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Hubbard model with two non-hybridizing bands: 

H = Yl ~ Amahma + A ^(n i2(T - h ila ) (2-1) 

<i,j>ma ia 

+ U ^ ^im^ n iml + nn a ni2a + U" ^ ^ila^a , 

im ia ia 

where c\ ma (cj ma ) creates (annihilates) an electron with spin cr(=j,J,) and orbital 
index m(= 1,2) at the ith site. The interorbital interactions U',U" are described by 
combinations of the on-site Coulomb interaction U and the local Hund's coupling 
integral J. We chose the parametrization U'=U — 2J and U"=U — 3J, proven to 
be suitable for t2g substatesP^J Crystal-field effects are taken into account via the 
energy splitting 2A between the orbitals. Spin-Flip and pair-hopping terms were 
neglected in (|2-1|) . For the numerical applications we choose the lattice to be the 
infinite-connectivity Bethe lattice giving rise to a semi-elliptical density of states 
(DOS) of bandwidth W m for band m. 

This model is treated within dynamical mean field theory (DMFT): the lat- 
tice model is mapped onto a local-impurity problem in an effective self-consistent 
bath. The self-consistency condition relates the local lattice Green's function to the 
impurity problem in the usual manner.^ The impurity problem is solved using 
the quantum Monte-Carlo (QMC) Hirsch-Fye algorithm^ with up to 128 slices in 
imaginary time r and 10 6 sweeps. 

In the following, n denotes the total filling, and n\,ri2 the respective subband 
filling. Figure 1 shows results for the two-band model with different bandwidths 
(ratio b=W2/W\=l/2) in the absence of crystal-field splitting. The overall filling of 
one electron (n=l) is kept constant during the calculation, but the partial fillings 
change from their non-interacting values of ni~2/3 and ri2~l/3. Figure la visualizes 
the non-interacting densities of states (DOS), while Fig. lb displays the partial 
occupations of the two bands as a function of the Coulomb interaction U for different 
values of U /J. The overall tendency to equalization of the orbital occupancies with 
increasing Coulomb interaction and in particular close to the insulating regime at 
large U is a consequence of the bandwidths differences being less important the 
stronger the renormalizations become. We note that the quasiparticle residues Z m , 
as extracted from the linear regime of the imaginary part of the self-energies on the 
Matsubara axis, are of similar magnitude (even if, interestingly the broad band is 
slightly stronger correlated, due to its bigger occupation). The influence of J which 
reinforces the orbital compensation is evident: for decreasing U/J ratios the orbital 
compensation effect is substantially more pronounced. 

For small J - the effect is striking at J=0 but also present for selected U/J - a 
small-C/ regime is observed where increasing mutual electronic interactions first lead 
to further orbital polarization. This should stem from the fact that by increasingly 
populating the broader band the electrons can gain kinetic energy and thus may 
compensate the increased potential energy. To quantify this subtle balance more 
rigorously we have considered the above model in a static approximation which allows 
for an exact solution. We replace the operators by their mean values h m — > {h m )=n m 
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Fig. 1. Two-band model calculations with W2=0-5W\, n=l and A=0. (a) Non-interacting DOS. 
(b) Band fillings for different U / J series in the metallic regime, (c) Band-resolved QP residue 
for different U/J series. Dark (red) lines belong to properties of band 1(2). 
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Fig. 2. Two-band model calculations in the incommensurate case with 1^2=0. 5Wi, n=0.5 and A=0. 
(a) Non-interacting DOS. (b) Band fillings for J=0 and (7/J=4 series in the metallic regime, 
(c) Band-resolved QP residue Z for both series. Dark (red) lines belong to properties of band 
1(2). 



and use a simple Friedel representation for the kinetic energy: £' kin ~n(2 — n). From 
a minimization of the sum of the kinetic and on-site potential energy one obtains 
the occupations 

{l + b(n-l))W 1 -nU{l-\) + 2A 

m = , n-j = n — ni . (2-2 

(1 + b)W x - 2U(1 - I ) ' 1 ' 

with n=n\ + n2 and a=U/J. Choosing b<l and refering to mo as the occupation in 
the non-interacting case, the condition for ni>ni j o reads: 

1. a>5 for 4>_i( n -2)(l-&) 

W-i 

2. a<5 for A < — -{n - 2)(1 - b) 
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Fig. 3. Two-band model calculations with Wi=W\, n=l and varying crystal field A. (a) Non- 
interacting DOS. (b) Band fillings for different U, J settings in the metallic regime, (c) Band- 
resolved QP residue Z. 

Independent of bandwidth, filling and crystal-field splitting, the system will thus 
polarize for U/J>5 and compensate for U/J<5, with respect to the non-interacting 
case. This limiting value can already be obtained while neglecting the kinetic- 
energy terms. It follows, e.g., from a Hartree-Fock treatment of the interaction 
term aloneEHEHID I n fact, it can be shown that in the static limit the polariza- 
tion/compensation regimes are in most cases independent of the specific form of the 
band energy. In the small-?/ regime where fluctuations are weak and the bandwidth 
renormalizations small, the static limit indeed describes the charge flow correctly (see 
Fig.l). However, for larger U dynamical effects and the changes in kinetic energy 
due to the bandwidth renormalizations become important: in this regime, the system 
tends to level out orbital occupation differences and eq. (|2-2j) no longer applies. 

The situation is different for non-integer total filling n, (see Fig. |2J). Though for 
small U the behavior is again qualitatively correctly described by (|2-2|) . the system 
always tends to polarize in the large-?/ limit. Indeed, in the doped case there is no 
Mott transition and fluctuations are weaker, even for large U (Fig. [2b)- Interestingly, 
there appears to be also a change in the degree of correlation for the different bands 
when turning on J. 

Finally, we present results from calculations for a two-band model with finite 
crystal-field splitting and equal bandwidths Wi=W2 in the commensurate case with 
n=l. Increasing Coulomb interactions drastically increase the orbital polarization, 
the stronger the larger the crystal field A. However, a large Hund's coupling J works 
against this polarization, as seen for U / J=4, which brings the orbital occupations 
back to their non-interacting values. As seen in Fig. |3b, the large differences in the 
occupations also give rise to substantially different quasiparticle residues Z m . 

From these model calculations we conclude that in general Coulomb interac- 
tions in non-degenerate multiband systems can induce substantial charge transfers 
between single subbands. Hund's rule coupling generally leads to a compensation 
of orbital occupation differences, whereas crystal- field splittings tend to polarize the 
system. The net charge flow depends on the fine balance of these effects, renormal- 
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ized by the Coulomb interactions and influenced by the anisotropics of the system. 
This explains, why indeed in some systems such as the series of 3d 1 perovskite com- 
pounds studied in Ref. [7J) correlations were found to orbitally polarize the system 
whereas in others (e.g. S^RuOj^'^) compensation effects were observed. Below 
we discuss two further recent examples, in which the question of orbital polariza- 
tion or compensation effects acquires deeper importance by the fact, that the charge 
transfers induce substantial Fermi surface modifications: the first one is the vana- 
dium sulfide compound BaVSs where the charge compensation effect was recently 
proposecP to be at the origin of a Fermi-surface modification that could be cru- 
cial for explaining the low-temperature properties of this compound. The second 
example is the thermoelectric cobaltate system Na x Co02, where very recent angle- 
resolved photoemission (ARPES) experimental dat^^'^ seem to hint towards a 
charge polarization effect, whereas there are conflicting results concerning this issue 
from theoretical investigations of this problem.^ 'E^JJnJ 

§3. The metal-insulator transition of BaVSs 

3.1. Brief overview of the physical properties of BaVSs 

BaVSs is well known as a long-standing challenging compound, exhibiting un- 
usual electrical and magnetic properties J21J2J JUDJH) At room temperature, 
this vanadium sulfide crystallizes in a hexagonal (PQ^/mmc) structurepS' whereby 
straight chains of face sharing YSq octahedra are directed along the c axis. Remark- 
ably, BaVS3 undergoes three distinct continuous phase transition with decreasing 
temperature. At T s ~240 K the crystal structure transforms into an orthorhombic 
(Cmc2i) structure (see Fig. |1J), 26! thereby creating an anisotropy in the o6-plane, 
i.e., perpendicular to the chain direction, and a zigzag distortion of the VS3 chains 
in the ftoplane. Additionally, the sign of the Hall coefficient changes from negative 
to positive.^ The orthorhombic as well as the hexagonal unit cell both contain two 
formula units of BaVS3, whereby the relevant V atoms are equivalent by symmetry. 
The system displays a MIT at T MIT ~70 K into a paramagnetic phase, hence BaVS3 is 
an example where a MIT occurs without intervening magnetic ordering. Apparently, 
a magnetic transition finally takes place at T x ~30 K, where an incommensurate an- 
tiferromagnetic order seems to be established!^ 

The paramagnetic conducting phase behaves as metallic (dp/dT>0) above ~150 
K, from where the resistivity^ (see Fig. |SJ) and the Hall coefficienlP^ ( se e Fig. EJ) 
increase up to the MIT. From diffuse scattering experiments,^ correlated struc- 
tural fluctuations along the direction of the VS3 chains were deduced. These fluc- 
tuations were interpreted as a precursor for a commensurate charge density wave 
(CDW) instability associated with a possible Peierls mechanism that triggers the 
MIT. The obtained wave vector for the structural instability in the orthorhombic 
cell Qmit— (1 5 0, \)o is identical with the one derived by Inami et a/. 29 -' which indeed 
detected a doubling of the unit cell below T MIT . In a recent x-ray study Fagot et 
were able to characterize the crystal structure below 70 K. The unit cell with a mon- 
oclinic distortion (space group Cm) now contains four formula units, whereby there 




appears to be a trimerization of the V atoms along the c axis, corresponding to a 
dominant 2fcp distortion.^ It should be kept in mind however that the conduction 
anisotropy within the system is not strongly pronounced (<7 c /<7 a ~3-4)p^ making 
pure ID interpretations questionable. Furthermore, the "bad metal" regime above 
the MIT with the resistivity minium at ~150 K is still not quite understood. In 
ultraviolet photoelectron-spectroscopic 33 ' and photoemission2=l*23 studies no Fermi 
edge was observed in the temperature regime around T M it This raised some specu- 
lation about a possible realization of a Luttinger electron liquid in this compound, 
which however seems unrealistic due to the intriguing interplay of 3D and ID fea- 
tures. The large increase of the Hall constanlP^ within the precursive regime seems 
to suggest that carriers are scattered into immobile states. In addition, an impor- 
tant feature of the phase above T MIT is the existence of local moments as revealed by 
the Curie- Weiss form of the magnetic susceptibility.^ The effective moment cor- 
responds approximately to one localized spin-1/2 per two V sites. Since the formal 
valence is V 4+ , corresponding to one electron in the 3d-shell, this can be interpreted 
as the effective localization of half of the electrons. At T MIT , the susceptibility rapidly 
drops, and the electronic entropy is strongly suppressed.^ 

The electronic structure close to the Fermi level in the orthorhombic (Cmc2\) 
phase is dominated by low-lying V(3d)-t2 g states. Due to symmetry, these states are 
split into one A\ g and two E g states per V atom. First-principles DFT-LDA cal- 
culations for Cmc2i-BaVSj^ do yield a V(3d)-S(3p) hybridization which is strong 
enough to account for the weak anisotropy of the transport properties. No band-gap 
opening has been reached within LDA. Instead, very narrow E g bands right at the 
Fermi level, and a nearly filled dispersive band with mainly A\ g character extending 
along the c* direction (c* is the reciprocal unit cell vector along the c axis of the 
system) have been found. This is consistent with a simple model proposed early 
on by Massenet et aiP^J However, the occupancy of the narrow E g bands found 
within LDA is too low to account for the observed local moment in the metallic 
phase. The nature of the CDW instability is also left unexplained by the LDA calcu- 
lations. Indeed, the calculated norm of the Fermi wave vector of the broad A\ g band 
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Fig. 5. Electrical resistivity (full line) and 
magnetic susceptibility (circles) of BaVS3 
(from Ref.H2l). 
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Fig. 6. Hall coefficient Rh versus tempera- 
ture. The large increase close to Tmit in- 
dicates the carrier loss (from Ref. \'2'Sl ). 



is found to be 2fc^P A ~0.94c*p^' while the observed wave vector of the instability 
is q MIT =0.5c*P^ Therefore, the simple picture of a CDW at q MIT =2kp associated 
only with the A± g band is untenable within LDA. It is likely that the E g states 
also participate in the instability, still the LDA band structure does not provide a 
Fermi-surface nesting that is in line with experimental findings. Hence, ab-initio 
calculations based on L(S)DA are not sufficient to explain the complex electronic 
structure of BaVS3. Static L(S)DA+U and/or GGA+U schemes appear also not to 
be adequate. Despite the fact that a band gap is obtained for Cmc2i-BaVS3 when 
enforcing magnetic order, this does not resemble the experimental scenario of an 
CDW mechanism underlying the transition from a Curie- Weiss-like metal into an 
paramagnetic insulator. 

In the following we will show that within a DMFT framework, using the LDA elec- 
tronic structure as a starting point, one may indeed reconcile theory with experi- 
mental findings. On the basis of this LDA+DMFT treatment we propose correlation 
effects in a multi-orbital context as an explanation for the discrepancies between band 
theory predictions and experiments, in line with the results of our model investiga- 
tions in 32 Specifically, we show that interorbital charge transfers occur which lower 
the occupancy of the A± g orbital in favor of the E^s. Thus leading to a Fermi surface 
modification by mainly shifting kp{A\g) towards lower values. From a calculation 
of the local susceptibilities, we demonstrate that local moments are formed in the 
metallic phase due to the low quasiparticle coherence scale induced by the strong 
correlations (in particular for the narrow E g bands). 

3.2. Calculational Approach 

We performed a realistic many-body investigation of the electronic structure of 
BaVS3 by means of the combination of the LDA with the DMFT. The ab-initio 
LDA calculations were performed with a pseudopotential code within a mixed-basis 
consisting of plane waves and localized functions.^ Norm-conserving pseudopoten- 
tials were used, and for the exchange-correlation functional the parametrization of 
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Fig. 7. Symmetry adapted V(3d) orbitals for BaVS3 in the Cmc2\ structure, (a) A\ g , (b) E g as 
well as (c) e g . In the orthorhombic Cmc2\ structure the following correspondance for the axes 
is realized: a «-> x, b <-» y, c <-» z. 



Perdew and Wang-- was utilized. In order to account for the given crystal sym- 
metry of BaVSs, a symmetry-adapted V(3d)-basis {4> m } (see FigJTJ) was obtained 
by diagonalizing the orbital density matrix umm"^ Skf> fkbi^kblM) {M' , where 
V'kfe stands for the pseudo crystal wave function for wave vector k and band b, and 
M,M' denote the cubic harmonics for 1=2. 

In the LDA+DMFT approach the on-site Green's function of a multi-orbital prob- 
lem reads as 

G a (iu> n ) = [(iw m + /i)l — H(k) — Xl (J (zu; ri )] _1 , (3-1) 

k 

with H(k) as the self-consistent determined LDA hamiltonian, provided in a localized 
basis, as well as the ^-independent self-energy ~E(iuj n ). In Eq. (|3- 1|) . a is the spin 
index, uj n ={2n + (3 (n = 0, ±1, ±2, . . .) are the fermionic Matsubara frequencies 
for the inverse temperature /?=T _1 , and fx is the chemical potential. For the special 
matter in hand, we restricted the calculations to the paramagnetic case, because no 
magnetic order is expected in the MIT regime. Furthermore, we approximated the 
full LDA hamiltonian by its downfolded version onto the {^4i 9 ,.E 9 } basis, which is 
sufficient to capture the essential physics close to the MIT. This downfolding to an 
effective 3-band model was performed 'empirically', lead by straightforward physical 
arguments, on the LDA DOSP The V(3ci)-Green's function matrix for the effective 
3-band model takes diagonal shape in the LDA limit. We kept this shape also within 
the DMFT framework, since off-diagonal self-energy terms should be small due to 
symmetry. Thus one may replace the /c-sum in eq. ()3-l|) by the integral over the 
partial DOS D m (e), so that the orbital-resolved Green's functions are written as 

f deD { m DA) (e) 
G m (iuj n )= / ■ — — — - , m = Aig,E g i,Eg 2 . (3-2) 

J lUJ n + fJ, ~ E - £ m {lUn) 



Note that by treating only correlated states in this realistic DMFT formalism, 
double-counting terms originating from correlations already included in the LDA are 
absorbed by the chemical potential [i. The on-site vertex was again parametrized as 
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U^ m =U, U^,,=U-2J and U^X^=U-3J, with U the on-site Coulomb repulsion 
and J the local Hund's rule coupling. Recall that J does not only describe the spin 
exchange energy, but also the reduction of U for electrons in different orbitals. The 
local impurity problem within the DMFT formalism was again solved using QMC 
with similar convergence parameters as in the model investigations. 

3.3. The electronic structure within LDA 

Figure El shows the computed band structure for the Cmc2\ crystal data at 
T=100 Kp§) The high-symmetry points T-C-Y define a triangle in the k z =0 plane, 
whereas Z-E-T is the analog shifted triangle in the A; 2 =0.5c* plane. The T-Z line 
corresponds to the propagation along the chain direction in BaVS3. Also drawn 
are the so-called fatbands for the {Ai g ,E g } orbitals. The width of these fatbands is 
proportional to the amount of the orbital character of a given band in the symmetry- 
adapted {4> m } basis. Being directed along the chain direction between pairs of V 
atoms, the A\ g orbital has mainly d z 2 character. In contrast, the E g states, linear 
combinations of d yz , d x %_ y i and d z i (E g i) as well as d xy and d xz (Eg?), only weakly 
hybridize with their surrounding. The orbitals of the remaining e g manifold point 
mainly towards the sulfur atoms, which results in a large energy splitting, leading to 
a smaller (larger) contribution to the occupied(unoccupied) states well below(above) 
ep- Therefore, the e g states do not have a major influence on the essential physics 
around the MIT. In Fig. |S] one can clearly identify that the very narrow bands at 
the Fermi level are associated with the two E g orbitals. Along T-Z starting at 
around -1 eV, one observes a nearly filled dispersive band with mainly A± g character 
that crosses the Fermi level close to the edge of the Brillouin zone (BZ). The 2kp 
value for this band amounts to 0.94c*, hence is incompatible with the experimentally 
observed 2kp=0. 5c* instability. In addition, the filling of the E g bands is rather low 
and cannot account for the observed magnitude of the local moment. Note however 




Fig. 8. LDA band structure for Cmc2i-BaVS3 
along high-symmetry lines in the 1. Bril- 
louin Zone. Also shown are the fatbands 
(see text) for the A\ g (blue), E g i (red) and 
E g 2 (dashed-green) orbital of the V atoms. 
The wiggling of the fatbands in the Z-E-T 
plane is due to the high degeneracy of the 
bands. 




E-E p (eV) 



Fig. 9. Partial LDA-DOS of the V(3d) states 
in the symmetry-adapted basis. For the 
computation of this partial DOS, a sphere 
of radius 2.18 a.u., i.e., half of the nearest- 
neighbor V-S distance, was introduced. 
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the small E g i electron pocket at the r point. It is absent in the band structure 
of the hexagonal phase and may be related to the hole-like conduction below T$ as 
seen in Hall measurements The relevant /c-summed electronic structure can be 
examined in Fig. |§] via the TO -resolved LDA-DOS per atom of the V(3d) states. 
This plot of the partial DOS is characterized around the energy intervall [-1,2] eV 
by a rather broad A\ g band and two strikingly narrow E g bands right at Ef, in 
accordance with the shown band structure. The small crystal-field splittings read: 
A(A lg ,E g )=e Alg - ££ 9 ~0.1 eV and A{E gl ,E g2 ) = e Egl - e Eg2 ~-0.05 eV. 

3.4. Orbital repopulations and Fermi surface deformations 

We will show that the intimate interplay between the Hubbard parameters U and 
J as well as the given LDA band dispersions can eventually result in a substantial 
charge transfer from the broader A± g band to the narrower E g bands. Additionally, 
accompanied by this charge transfer is a shift of the A\ g band in k-space, providing 
the possibility for the onset of the experimentally detected CDW instability. 
Due to the lack of information about the value of U in BaVS3 from experiment 
(e.g., photoemission) or theory (constrained LDA methods tend to underestimate the 
screening for metals), we undertook our investigation by carefully scanning through 
the physically meaningful U-J parameter space for this compound. Since the in- 
terplay of U and J appears to be a central issue, we chose to fix the ratio U /J. 
By varying U we studied two different series: U/J=7 and U/J=A. Recall that in 
the 2-band model we encountered in £JH the ratio U / J=5 separated the polariza- 
tion/compensation regimes in the static limit. One would expect to find similar 
regimes also in the present 3-band model for BaVS3. 

The orbital occupancies in our effective 3-band model, at the LDA level (i.e., for 
U=0) read: n{A lg )=0.7U, n(E gl )=0.207 and n(E g2 )=0.081. In Fig. [[UK we plotted 
the respective band fillings in the metallic regime for increasing U within the two Uf J 
series. The main effect apparent on Fig. HOh . is that moderate correlations indeed 
tend to bring the occupancies of each orbital closer to one another, i.e., to decrease 
the population of the "extended" A\ g orbital and to increase the occupancy of the 
E g orbitals. For strong correlations, values close to n{^A\ g )~n{E g \)+n{E g 2)— 0.5 are 
obtained in the DMFT calculation, corresponding to a half-filled band. The differ- 
ence in the qualitative [/-dependence of the band fillings between the two U/J series 
resembles those of the model calculations. For U/J=4 there is an immediate orbital 
charge transfer for Uj^O, whereas for U/J=7 correlations first have to overcome the 
polarizing tendencies for small U. Hence the basic mechanisms found in the model 
calculations appear to be an essential point in understanding the physics of BaVS3. 
Note that observations on the importance of J and of band-narrowing effects for 
the band filling in realistic systems have also been pointed out in the context of 
ruthenatesP^ For C/~4 eV the system reaches the Mott transition within U / J=7, 
whereas this transition is shifted to much larger U values for U/J=4. In the present 
context, we did not investigate the Mott insulating state as this one may not be rele- 
vant for realistic BaVS3. Recall that in nature the orthorhombic (Cmc2i) structure 
transforms at T MIT via an CDW instability into a monoclinic (Cm) structure with 
4 inequivalent vanadium atoms in the unit cell. Hence an understanding of the true 
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Fig. 10. (a) Band fillings at /3=(fc B r) _1 =15 
eV" 1 (T=774 K) for the effective bands 
within LDA+DMFT. (b) Corresponding 
shift of the Fermi level for the A\ g band 
(note that '=0). Filled symbols: 

U/J=7, open symbols: U/J—4. 
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Fig. 11. LDA+DMFT spectral data for (7=3.5 
eV, U I J=4. (a) V(3d) low-energy quasi- 
particle bands along F-Z in LDA (dashed 
lines) and LDA+DMFT (solid lines) for 
T=332 K. (b,c) integrated spectral func- 
tion p{uj) for a single formula unit of BaVS3 
at T=1160 K (b) and T=332 K (c). 



insulating state cannot be obtained by driving the system in the Mott-insulating 
regime. For the strongly correlated metal close to T MIT we would estimate C7~3.5 
eV in order to explain the substantial filling of the E g states revealed from the local 
moment examinations. 

As already noted, the correlation-induced orbital charge transfers not only modify the 
respective band fillings, but may lead also to relevant changes of the Fermi surface. 
Such an effect is anticipated in the case of BaVS3 due to the insufficiency of DFT in 
LDA to explain the observed CDW instability. Indeed our calculations reveal that the 
depletion of the A\ g band is accompanied by a reduction of the corresponding Fermi 
wave vector along the F-Z direction. While a full determination of the quasipar- 
ticle (QP) band structure in the interacting system requires a determination of the 
real-frequency self-energy, we can extract the low-energy expansion of this quantity 
from our QMC calculation in the form: KeU rn (uj+iO + )c^ReU m (0)+uj(l — 1/Z m ) + - ■ ■ , 
with Z m the QP residue associated with each orbital. The poles of the Green's func- 
tion determine the QP dispersion relation: detfwk — Z[H^ DA + Rei7(0) — /i]]=0, with 
[i the chemical potential. Focusing first on the A\ g sheet of the Fermi surface, within 
our diagonal formulation the location of the Fermi wave vector in the interacting sys- 
tem is determined by: e^ A (ki?)= // — KeHA lg (0). This quantity therefore yields the 
energy shift of the A\ g band at the Fermi surface crossing, as compared to LDA. It 
is depicted in Fig. HUb as a function of U. In Fig. [TTk . we display the QP bands that 
cross the Fermi level along F-Z in a narrow energy range around ef for U =3.5 eV in 
the U I J=4 series. The QP bands are obtained by performing a perturbative expan- 
sion of the pole equation above, which yields: oj^= J2 m C^ nk Z m [e^ A +KeE m (0) — //] 
with C^KV-kbl^m)! 2 the LDA orbital weight. From these two figures, it is clear 
that k,F{Aig) is reduced in comparison to the LDA value, in line with the global 
charge transfer from A\ g to E g . This opens new possibilities for the CDW instabil- 
ity, in particular for the nesting wave vector. 
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Fig. 12. Temperature-dependent local spin susceptibilities for [7=3.5 eV, according to the normal- 
ization x' 1 ° c ' 1= /o 3 dT{S z {0)Sz{T)) , where S z denotes the z-component of the spin operator. 



The enhanced population of the narrow E g bands, as well as the correlation-induced 
reduction of its bandwidth (see Fig. II lb ) provide an explanation for the local mo- 
ments observed in the metallic phase. To support this, we have calculated (Fig. 112)) 

the local susceptibility associated with each orbital Xm^ = J2q^ e [Xm(<l, lo=0)]. For 
both values of U / J, the susceptibility of the A\ g band saturates to a Pauli-like value 
at low temperatures. In contrast, x^ loc ^ °f the- E g orbitals strongly increases as the 
temperature is lowered (except for the low-filled E g 2 orbital at U/J=A). This is 
because the coherence temperature below which quasiparticles form is much lower 
for the E g orbitals than for the A\ g orbital. Accordingly, our calculation of the in- 
tegrated spectral functions (Fig. ITTb.c) reveals a strong T-dependence of the E g QP 
peak. Some differences between the two series are clear from Fig. 1121 For U /J=7, 
the system is already very close to the Mott transition when choosing [7=3.5 eV. 
Thus the A\ g electrons also act as local moments over part of the temperature range, 
while for U/J=4 and same U, the T-dependence of the total local susceptibility is 
almost entirely due to the E g \ electrons. Which of the two situations is closest to 
the physics of BaVS3 does require further investigations, albeit some experimental 
indications point at the second possibility.® 

§4. The Fermi surface of Na^CoC^ 

Since the surprising discovery of superconductivity in NazCoCvyH^C® for 
x~0.35, the sodium cobaltate system has made its way to one of the hottest topics in 
current research on strongly correlated materials. The layered C0O2 oxides, held to- 
gether by Na layers in between, resemble the famous cuprates, although the Co atoms 
form a triangular lattice with the oxygen atoms in a pyramidal configuration above 
and below the Co plane. Besides the superconducting phase the unhydrated sodium 
cobaltates exhibit additional rich physics, such as strong electric thermopower,® 
charge order,® different magnetic order,® as well as complex ordering of the Na 
atoms.® Bounded by a Mott-Hubbard insulator for x=0, i.e., C0O2, and by the 
band insulator NaCo02 for x=l, the sodium cobaltate phase diagram® is roughly 
separated in a Pauli-like metallic region for x<0.5 and a Curie- Weiss-like metallic 
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region for x>0.5. In another difference to the cuprates, cobaltate is a multi-orbital 
3e? 5+x system. Due to the crystal-field splitting the Co atoms are in a low-spin state, 
hence the d electrons occupy the lower t2 g band complex (i.e., there are 1 — x holes 
in the t2 9 shell). It follows that the Curie- Weiss behavior stems from Co 4+ local 
moments, whereas these moments seem to disappear below x=0.5. Interestingly, 
magnetic phases are experimentally verified for x>0.75. 46 ' 

From this short introduction of Na x Co02 it is already obvious that the x=0.5 
composition plays a subtle role in this system. Although the characteristics of the 
phase at x =0.5 are still not completely resolved, certain properties seem to be 
established. In contrast to most other doping levels x, Nao.sCoC^ exhibits an ordered 
superstructure for the combined Na/Co02 system^ over a wide temperature range. 
Most importantly, Nao.sCoC^ shows a MIT at ~53 K. If this transition is truly 
associated with the onset of charge ordering is still a matter of debate, although 
many measurements point in this directionP^J^ There appear to be^ additional 
magnetic and/or structural transitions at 20 K and 87 K. 

Early electronic structure calculations in the L(S)DA framework for Nao.sCoC^ 
have been performed by Singh. 49 ' Since the the ordered superstructure was identified 
only recently, in these DFT calculations a disordered arrangement of the sodium 
atoms was assumed (treated within the virtual crystal approximation). The Fermi 
surface revealed in this investigation exhibits a large cylindrical hole sheet with 
dominant A\ g character around the r~A line in the Brillouin zone of the hexagonal 
unit cell, as well as small hole pockets in the F-K and A-H directions. This topology 
of the Fermi surface appears to be rather generic for a wide range of x within the 
LDA framework. It has been sugg estecPJ that strong nesting between the hole 
pockets leads to large fluctuations in the spin channel, having important influence 
on the character of the superconducting state at x~0.35. However the existence of 
these hole pockets has not been confirmed in angle-resolved photoemission (ARPES) 
experiments JaDflsuCD On the contrary, in one of these ARPES measurements- 1 -*' the 
"pocket" bands were found well below the Fermi level for large doping range x. 

It appears obvious from the latter discussion that Fermi-surface modifications 
due to strong correlation effects may play also an important role for sodium cobal- 
tates. Hence the last sections are devoted to a brief analysis of the problem of finding 
the correct generic Fermi surface for Na x Co02- 

4.1. Electronic structure within LDA 

In order to have a well-defined starting point, we tried to find theoretical agree- 
ment with the experimentally detected ordered superstructure for Nao.sCo02. Thus 
we performed total energy calculations with the mixed-basis pseudopotential (MBPP) 
code for different arrangements of sodium atoms in the interlayers in between the 
C0O2 planes. As a starting point for the Na x Co02 structural model one usually 
utilizes the hexagonal (P6322 or PQ^/mmc) symmetryP^^^J^ In all our studies 
the lattice parameters were confined to the values obtained by Jansen and Hoppej^ 
i.e., a#=2.84Aand c#=10.8lA. Indeed we obtained the lowest LDA total energy 
for the experimentally suggested superstructure with global orthorhombic (Pmmn) 
symmetry.^ We relaxed the atomic positions within the MBPP and the corre- 
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sponding superlattice is depicted in Fig. 1131 The unit cell of this superstructure 
with dimension ao=2aH ,bo = aiiV3 and co=ch contains 28 atoms, whereby there 
are 2 inequivalent Na and Co sites, respectively, as well as 3 inequivalent O sites. 
The two classes of inequivalent Na and Co sites are related as the Na ordering is such 
that in each sodium interlayer there are atoms sitting on top of an Co site (Nal/Col) 
and such atoms that are placed on top of an oxygen site (Na2) (see. Fig. IT3b). The 
Co atoms without an Na atom on top form the Co2 class. As reported by Huang et 
aJ.j^ this ordered superstructure should be stable over a wide temperature range, 
with possible structural changes below 100 K. 

We determined the LDA electronic structure for the Pmmn symmetry with 
the MBPP code (a similar study was published recently^). The resulting band 
structure and the summetry-adapted 3d-DOS for the two Co classes are shown in 
Fig. El an d Fig. 1151 The band structure around the Fermi level consists of 3 iso- 
lated blocks. Roughly speaking, the block deep in energy is formed by the 0(2p) 
states and the one high in energy stems from the Co(e 5 ) states. The most inter- 
esting block of bands close to the Fermi level is mainly dominated by the Co(i2g) 
multiplet and has only minor 0(2p) weight due to the remaining Co(t2 g )-0(2p) hy- 
bridization. Thereby, the Ai g subband complex is characterized by an substantial 
bonding-antibonding splitting, in between which the E g subband complex is mainly 
located. A crystal-field splitting of A=EA lg ~ £ e 3 ~-0.1 eV between these two band 
complexes can be derived. The crystal-field splitting between the different E g states 
appears to be negligible. From Fig. El no striking difference between the electronic 
states within the named t2 g multiplet for the two Co classes, i.e., Col and Co2, can 
be extracted. For the respective LDA orbital occupations the following inequalities 
hold: n^" 1 ^ >n^° and nE gll E g2 > n A lg - Hence the two Co classes are discriminated 




Fig. 13. Nao.5Co02 in the orthorhombic (Pmmn) structure, (a) 3D view, (b) view along the c 
axis. 
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Fig. 14. LDA band structure for Nao.5Co02 
in the Pmmn crystal structure along high- 
symmetry lines in the 1. BZ of the or- 
thorhombic lattice. 
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Fig. 15. Partial LDA-DOS of the Co(3d) 
states in the symmetry-adapted basis for 
Col and Co2 atoms. 



in LDA by a minor charge disproportionation (~0.06 electrons), and the results sup- 
port the suggestion™ that the Col atoms should tend to the lower oxidation state, 
i.e., Co 3+ rather than Co 4+ . However in this work we do not address the question of 
charge ordering. Finally, the hole occupancy is enlarged in the A\ g orbital for both 
Co classes. 

For the discussion of the LDA Fermi surface we plotted in Fig. ^] also the band 
structure close to the Fermi level for the larger (a=a#) hexagonal BZ (with symme- 
try lines in direct comparison to Singh's worlP^), since the generic Na x Co02 system 
is most often associated with the hexagonal symmetry stemming from the C0O2 sub- 
lattice. In this representation higher BZs of the orthorhombic structure are touched, 
thus the band stucture exhibits some folding. We also provide the fatbands for the 
respective t2 g orbitals in Fig. 1171 By comparing the band structure of the ordered 
superstructure in Fig. ED with the one assuming disordered sodium interlayers in a 
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smaller unit cell from Ref. 49), some differences may be identified. There are bands 
in the A-T direction right at the Fermi level which result in an additional Fermi sur- 
face sheet for the orthorhombic structure. ARPES intensity close to the r point has 
been detected recently but below the Fermi levelP^ In the same ARPES study the 
famous pocket bands were found below the Fermi level, whereas in our LDA calcula- 
tion they form, of course, still hole Fermi sheets. It is instructive to study the main 
orbital character of the relevant bands within the t2 g block. From Fig. I17| the hole 
pockets are nearly exclusively stemming from the E g orbitals of Co. The same holds 
for the bands along A-T at the Fermi level. As a general subtle difference between 
Col and Co2, one observes that in the case of Col the E g bands have either E g \ or 
E g 2 character, whereas for Co2 these bands have mixed E g x t 2 character. Since the 
large hole Fermi sheet has dominant Ai g character ARPES measurements suggest 
a strong orbital polarization for Nao.5Co02, namely with nearly entirely filled E g 
bands and A\ g hole bands with a filling close to 0.5 holes per Co atom. 

4.2. Pockets or no pockets? 

Since the i2 9 -manifold bandwidth W is very narrow, Na^OoC^ is expected to 
be in the regime of strong correlations W<tiU. This is qualitatively consistent with 
the substantial mass renormalization and high thermopower. A value of t/~4 eV 
has been estimated.^ Hence, considerations from DFT-LDA may not provide the 
correct picture of the electronic structure, even for the Fermi surface itself. In a 
recent LDA+U calculatiorP^ the pocket bands were shifted below the Fermi energy 
for various dopings x (corresponding to increased orbital polarization), in agreement 
with the ARPES measurements. We have also performed (unpublished) LDA+U 
calculations of orthorhombic Nao.sCoC^ and confirm this result. However, the static 
LDA+U approach is certainly questionable for the metallic regime of Nao.5Co02- 
Moreover, the precise shape of the Fermi surface in LDA+U depends on the specific 
choice of the double-counting correction terms. Recently, an LDA+DMFT calcula- 
tion^ obtained the opposite effect for realistic values of J, namely a further stabi- 
lization of the hole pockets compared to LDA. In contrast, another study based on 
the Gutzwiller approximation reported that the hole pockets were shifted below the 
Fermi levelf^ but this is perhaps not surprising since U -^oo and J=0 were used in 
this calculation, hence favoring orbital polarization. Note that the sodium cobaltates 
have non-integer total hole filling (1 — x)e[0,l] per Co atom. As explored in 32 the 
direction of the charge flow leading to orbital compensation or polarization is the 
result of a delicate balance between different parameters, namely: the total filling, 
the Hund's rule coupling, the crystal-field splitting, and the bandwidth of each band 
(taking into account the renormalization of the latter two by strong correlations). 
Further studies are therefore necessary to fully settle the issue for Naa;Co02. To 
conclude, let us emphasize that the following aspects, neglected in previous studies, 
might be important in order to reach this goal: 

• Non-diagonal terms S mm i in the self-energy, renormalizing the interorbital hy- 
bridizations and therefore influencing charge flows. 

• Spin-flip and pair-hopping (non-Ising) terms in the Coulomb vertex. 

• Importantly, an orbital-independent double-counting correction within the \>2g 
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multiplet has usually been assumed when including dynamical correlations. 
While orbital-dependence of this correction is expected to be rather small 
and hence should not affect significantly our conclusions for the larger charge- 
transfer effects in the case of BaVS3, it may play a substantial role in the more 
delicate case of Na^OoO^- 
• Finally, non-local effects might be important, such as those induced by a strong 
interatomic Coulomb interaction. The latter have been proposed^ to be im- 
portant for the cobaltates, as also suggested by the observed tendency to charge 
ordering. 
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